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INTRODUCTION 

The purpose of this study was to select monthly weatner variables tnat coulo 
oe used to estimate yields of the Argentina corn crop. Tne corn-growing area 
in Argentina is indicated in Figure 1. Most of Argentina's corn is grown in a 
concentrated area in central Argentina where the climate is predominantly humid 
subtropical. The western edge of the soybean area in central Cordoba is semi- 
arid with warm to hot summers. There drought and high temperatures can be a 
problem during the growing season which begins with planting in October and 
November and extends to harvest in April through June. 

METHOD 

Multiple regression analysis of yield with selected agroclimatic indices 
was used to derive a suitable model. The regression equation for the corn model 
derived is: 

Y *oc+ BjT + B 2 (ETMETHi) +83 ( TX, ) + E 

where 

A 

Y = Estimated yield 

CL = Constant 

Bj * Coefficients of the variables j * 1 - 3 

T = Trend 1965-1980 

A A _ 

ETMETHi * ET - ET for month i: where ET * K*PET and K * ET/PET. 

TX.j * Maximum temperature for month i, and 

E * Unexplained error. 

\ 

Trend between the years 1965 and 1980 was determined from a plot of yields as 

shown in Figure 2. ET (evapotranspi ration) minus ^T (climatically appropriate 

evapotranspi ration) was used as an index representing soil moisture for plant 

growth. A more complete definition is given in the Appendix. Large positive 
A 

ET - ET values suggest wet conditions. 
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CORK PRODUCTION 


CORN YIELDS 


Figure 2. Plots of Production and Yield Versus Year for 
Argentina Corn. 
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In developing the models, various procedures of the Statistical Analysis 
System fSAS Institute, Inc., 1979) were used. The procedures used and the 
operations per:o r med with eacn are sunmarizec in tne Appendix. Tne selectee 
model r.ac tne mgnest ft 2 ana included variables that were significant at the 
10 per cent level and were agronomical ly meaningful. 

DATA 

Crop data for Argentina from 1960 to 1979 was obtained from the Foreign 
Agricultural Service (Alan Vandagrith, personal communication, 1982). A data 
set was created with year of yield as year of planting. Since the growing 
season spans two numerical years, the crop data was expressed in terms of "year 
+ 1" (or harvested year) and the meteorological data corresponding to any yield 
included data for that year and "lagged" data for the previous year. 

The meteorological data was created using the general Argentina 
meteorological station file. This file was composed of data from several 
different sources, including the Monthly Climatic Data for the World and the 
Servicio Meteorological Nacional in Argentina (R.E. Jensen, C.M. Sakamoto, and 
S.E. Mummert; August, 1974). Stations inside the corn-producing area were 
averaged by province and weighted according to the percentage that their pro- 
vince contributed to the entire country's production. Figure 3 shows the corn- 
growing area in Argentina and associated stations; Table 1 lists the stations 
used and their weights. 

The years between 1961 and 1980 (year of harvest) were used in the model 
since these years contained the most complete meteorological data. 

PROCEDURES 

Weather variables were selected from correlations and plots with yield. 
These variables, along with TREND, were tried in regression equations. The 
limited number of data years severely limits the number of variables that 
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SI AT ION NAME WEIGHT 


BUENOS AIRES .50 

Pergamino 

Junin 

Nueve de Julio 

CORDOBA .27 

Rio Cuarto 
Bell Ville 

SANTA FE .23 

Rosario 

Casilda 

Las Delicias 

Parana 

Pergamino 




Table 1. Meteorological Stations Used 
for the Argentina Corn Model. 
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can be Included in a regression model by limiting the number of degrees of 
freedom. This influenced the choice of model considered "best" for soybeans. 
Tnree weather variables plus TREND were considered the optimum model with 19 
years of data. Tne tnree weather variables are: 


A 

ETMETH11 ET - ET for November 

TX12 Maximum temperature for Oecember 

TX1 Maximum temperature for January. 


The coefficient for the November variable was positive which reflects 
moisture requirements just after planting. The coefficients for maximum 
temperature for December and January were negative indicating that temperatures 
too high in the peak summer months during the tassel ling and silking stages lead 
to stress on the plant. There is a critical growth period for corn. The 
statistics of the selected model are summarized in Table 2. 

TEST RESULTS 

A jackknife test was run on the model. In this test a year was eliminated 
from the crop data and the model was used to predict that year's yield. This 
process was dun*, for each successive year beginning with 1961. The results were 
reasonable with the greatest difference (3.24 quintals) between predicted and 
actual yield in 1976. That year severe drought badly damaged the corn crop and 
harvest was hampered by rains. Test results are printed on Table 3 and plotted 
on Figure 4. 
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APPENDIX 


Definition of Variables 

The difference between E7 and ET is an index of tne amount of moisture 
available for plant growth. Soil moisture depletion is based on evapo- 


transplratlon (ET) estimates, determined as follows: 

(ET) n « [ { (PET) n - ( p) n ] + ( P) ] 

AVO n J 

where 


I 


t 


’ * 


(ET) n * "Actual" evapotranspi ration, 

(S) n .i * Available moisture at end of n-1 month, 

AWC * Maximum water holding capacity, 

(P) n = Precipitation for month n, and 

(PET) n * Potential evapotranspira^lon for month n. 

A 

ET - ET measures the difference between the actual evapotranspi rati on and the 
"climatically appropriate" evapotranspiration, giving an indication of soil 
moisture supply and demand. 


Statistical Analysis System Procedures Used 

PROC CORR Computes correlation coefficients between 

variables, including Pearson product-moment 
and weighted product-moment correlation. 


PROC PLOT 


Graphs one variable against another, producing 
a printer plot. 


PROC STEPWISE Provides five methods for stepwise regression. 

Stepwise is useful when selecting variables to 
be Included In a regression model from a 
collection of Independent variables. 


PROC STEPWISE FORWARD Begins by finding the one-variable model that 

produces the highest r 2. For each of the other 
independent variables, FORWARD calculates F- 
statistlcs reflecting the contribution to the 
model if the variable were to be Included. 


Anfe* ’ 
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PROC STEPWISE BACKWARD 


pro: stepwise stepwise 


PROC STEPWISE MAXR 


PROC PETM 
PROC Z INDEX 


Begins by calculating statistics for a model 
including all the independent variables. The 
variables are deleted from the model one by 
one until all tne remain r » variables produce 
F-statisci:s signifies'’: a: tie .lo level. 

Tne stepwise metnod is a modification of tne 
forward selection tecnnique, differing in 
that variables already in the model do not 
necessarily stay there. After a variable is 
added (as in the forward selection method) the 
stepwise method looks at all the variables 
already included in the model and deletes any 
variable that does not produce an F-statistic 
significant at the .10 level. Only after this 
check is made and the necessary deletions 
accomplished can another variable be added to 
the model. ► 

(Maximum r 2 improvement) Unlike the three 
techniques above, this method does not settle 
on a single method. Instead it looks for the 
"best" two-variable model, the "best" three 
variable model, and so forth. 

Uses latitude and mean monthly temperature to 
calculate Thornthwaite's potential evapo- 
transpi ration for each month. 

Uses monthly PET's, prt ipitation, SS 
(beginning moisture in surface layer), 

AWCS (available water capacity in surface 
layer, SU (beginning moisture in the 
underlying layer), and AWCU (available water 
capacity in the underlying layer) to calculate 
Palmer's sot) moisture budget, drought index 
Z, ET, and ET. 
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